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ABSTRACT 

The masses of supermassive black holes in broad-line active galactic nuclei (AGNs) can be measured through 
reverberation mapping, but this method currently cannot be applied to very large samples or to high-redshift 
AGNs. As a practical alternative, one can devise empirical scaling relations, based on the correlation between 
broad-line region size and AGN luminosity and the relation between black hole mass and bulge stellar velocity 
dispersion, to estimate the virial masses of black holes from single-epoch spectroscopy. We present a revised 
calibration of the black hole mass estimator for the commonly used H/3 emission line. Our new calibration takes 
into account the recent determination of the virial coefficient for pseudo and classical bulges. 


Subject headings: galaxies: active — galaxies: nuclei 
general 

1. INTRODUCTION 

Supermassive black holes (BHs) play a fundamental role in 
many aspects of contemporary extragalactic astronomy. How¬ 
ever, direct methods to measure BH masses currently are lim¬ 
ited to very nearby, largely inactive galaxies. Reverberation 
mapping (RM; Blandford & McKee 1982; Peterson 1993) of 
broad-line active galactic nuclei (AGNs) provides a means of 
measuring the size {R) of their broad-line region (BLR), which, 
in combination with the velocity widths (AV) of the broad 
emission lines, yields a virial estimate of the BH mass 

R(AVf 

Mbh(RM) = /^-^. (1) 

The virial coefficient / depends on the kinematics, geometry, 
and inclination of the BLR. In practice, it is set by normalizing 
the RM AGNs to the correlation between BH mass and bulge 
stellar velocity dispersion (Mbh-ct* relation; Eerrarese & Mer¬ 
ritt 2000; Gebhardt et al. 2000; see Kormendy & Ho 2013 for a 
review) established by local galaxies with direct BH mass mea¬ 
surements (e.g., Onken et al. 2004; Woo et al. 2010). De¬ 
spite its utility, RM requires time-consuming, long-term spec¬ 
troscopic monitoring, and, to date, only ^ 50, mostly nearby (z 
<0.1) AGNs have been studied in this manner (e.g., Peterson 
et al 2004; Barth et al. 2011; Grier et al. 2013). To access 
more luminous quasars and especially higher redshift systems, 
an alternative, much more expedient mass estimator can be de¬ 
vised from single-epoch spectroscopy by taking advantage of 
the empirical correlation between BLR size and continuum lu¬ 
minosity, RocL^ (Kaspi et al. 2000), to estimate R. Then, from 
a straightforward measure of L and AV from single-epoch or 
mean (time-averaged) spectra of the RM sample, one can solve 
for 

logMBH(RM) = logVP(mean)-l-fl, (2) 

where a is a zero point offset and we define the virial product 
as 


galaxies: Seyfert — quasars: emission lines — quasars: 


/ AV \ V ALa \''' 

The line width AV can be parameterized as the full-width at 
half maximum (EWHM) or line dispersion (uiine, second mo¬ 
ment) of the broad line profile. The above approach to calibrate 
“single-epoch” spectra has been developed by a number of au¬ 
thors (e.g., McLure & Dunlop 2001; McLure &. Jarvis 2002; 
Vestergaard 2002; Vestergaard & Peterson 2006; Wang et al. 
2009). Eor lower redshift sources (z -^0.15), the most com¬ 
monly used emission line is H/3, and the continuum luminosity 
is referenced to 5100 A. 

Some recent developments justify a reassessment of the virial 
formalism to estimate BH masses from single-epoch H/3 spec¬ 
troscopy. Kormendy &. Ho (2013) significantly updated the 
Mbh - cr* relation for inactive galaxies, highlighting, in particu¬ 
lar, the large and systematic differences between the relations 
for pseudo and classical bulges. This prompted Ho & Kim 
(2014) to calibrate the / factor separately for the two bulge 
types, using the latest sample of RM AGNs for which reliable 
bulge classifications could be performed. The changes are not 
negligible. Whereas previous studies obtained / Ri 4.2-5.5 
(Onken et al. 2004; Woo et al. 2010; Park et al. 2012; Grier et 
al. 2013), Ho & Kim find / = 6.3 ±1.5 for classical bulges and 
ellipticals and / = 3.2 ± 0.7 for pseudobulges. (Eor the present 
discussion, we refer to / calculated using AV = criine-) Apart 
from these revisions to the Mbh - cr* relation and the new cali¬ 
bration of the / factor for pseudo and classical bulges, the R-L 
relation itself has been updated by Bentz et al. (2013). 

2. CALIBRATION 

Our new calibration of Equation 2 uses the updated database 
of properties for the RM AGNs and the bulge type classifi¬ 
cations of their host galaxies given in Ho & Kim (2014), as 
summarized in Table 1. Specifically, we compute Mbh(RM) = 
/VP(CTiine), using / and VP(CTiine) derived from root-mean- 
square (rms) spectra. We adopt / = 6.3 ± 1.5 for classical bulges 
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Fig. 1.— Comparison between RM BH masses and virial products calculated from FWHM of H/3 extracted from mean spectra for (left) the entire sample, 
(middle) classical bulges, and (right) pseudobulges. The solid line denotes a slope of unity, while the dashed and dotted lines give the regression based on the 
FITEXY and BCES fits, respectively. 


and ellipticals and / = 3.2 ± 0.7 for pseudobulges. As an ap¬ 
proximation to quantities pertaining to single-epoch spectra, we 
use the values of ALa( 5100 A) and FWHM(H/3) derived from 
mean spectra, as listed in Table 2 of Ho & Kim (2014). Of 
the 43 RM AGNs that have bulge type classifications, 38 have 
available H/3 line widths extracted from mean spectra, among 
them 14 pseudo and 24 classical bulges. Among the empirical 
several criteria recommended by Kormendy & Ho (2013; see 
their Supplemental Material) to distinguish pseudobulges from 
classical bulges, we adopt, whenever possible, the most widely 
used condition; Sersic (1968) index n <2. However, the Ser- 
sic index of the bulge can be difficult to measure accurately in 
the presence of a bright nucleus. Under these circumstances. 
Ho & Kim (2014) use a surrogate, equally effective criterion, 
that the bulge-to-total light fraction should be ^1/3. Six ob¬ 
jects have more than one set of observations; in our analysis we 
treat these multiple data points as independent measurements. 
We have verified that combining the multiple measurements as 
a weighted average does not alter our final conclusions. 

Figure 1 confirms that Mbh(RM) correlates tightly and lin¬ 
early with VP(mean). A linear regression using the FITEXY 
estimator (Press et al. 1992, as modified by Tremaine et al. 
2002) yields a formal slope of 0.97 ±0.07 and an intrinsic scat¬ 
ter of 0.35 dex for the entire sample with both bulge types 
combined. The BCES algorithm of Akritas & Bershady (1996) 
gives a consistent slope of 1.00 ±0.10. The middle and right 
panels of the figure show the correlations for the two bulge 
types separately. Although the FITEXY fit seems to suggest 
that pseudobulges have a slope less than unity, this is not sup¬ 
ported by the BCES fit. Overall there is no evidence that the 
two bulge types behave differently. 

Adopting a slope of 7 = 0.533^o 033 for th® latest R-L relation 
(Bentz et al. 2013), the final mass scaling relation becomes 


logMBH(H/3) = log 


/FWHM(H/3y /ALa(5100) 
\1000kms“'/ \ 10‘*"^ergs“* 


0.533" 




(4) 

where a = 7.03 ± 0.02 for classical bulges and a = 6.62 ± 0.04 
for pseudobulges, with a corresponding intrinsic scatter of 0.32 


and 0.38 dex. For both bulge types combined, a = 6.91 ±0.02 
with an intrinsic scatter of 0.35 dex, essentially identical to the 
results of Vestergaard & Peterson (2006). It is remarkable that 
our new fit, which is based on a larger sample and significant 
updates to all of the RM and velocity dispersion data, turns out 
to be so similar to that of Vestergaard & Peterson published 
almost a decade ago. This may indicate that this type of cali¬ 
bration is still currently dominated by systematic effects (e.g., 
intrinsic differences in the factor /). 

3. CONCLUDING REMARKS 

We present a new calibration of the prescription for estimat¬ 
ing BH masses for broad-line AGNs and quasars using single¬ 
epoch spectra of the H/3 emission line. The primary difference 
between this and the previous work of Vestergaard & Peterson 
(2006) is that we account for the systematic difference in the 
virial coefficient / between pseudo and classical bulges (Ho 
& Kim 2014). We explicitly assume that AGNs hosted by 
pseudobulges follow a different Mbh - cr* relation than those 
hosted by classical bulges, as observed in inactive galaxies. In 
addition, our calibration sample of RM AGNs is significantly 
larger and more current than that used by Vestergaard & Peter¬ 
son, which was based on the sample of Peterson et al. (2004) 
from a decade ago. Incorporating also recent minor changes to 
the R-L relation, we find that the zero point of the H/3 virial 
mass formalism for classical bulges is a factor of 2.6 (0.41 dex) 
higher than that for pseudobulges. The commonly used H^ cal¬ 
ibration of Vestergaard & Peterson (2006) is very similar to ours 
for classical bulges, but their zero point is a factor of 2 higher 
than our zero point for pseudobulges. 

This important source of systematic uncertainty obviously 
should be eliminated to the extent possible. However, as dis¬ 
cussed in Ho & Kim (2014), in practice this will prove challeng¬ 
ing because of the difficulty of measuring accurate host galaxy 
parameters for luminous, and especially distant, AGNs. As a 
general rule of thumb, any system with Mbh ^10* Mq can be 
safely regarded as a classical bulge or elliptical galaxy. This can 
be seen from the distribution of BH masses in inactive galaxies 
(Ho & Kim 2014, Fig. 2), as well as for the RM-mapped AGNs 
in this paper (right panel of Fig. 1). At the other extreme, BHs 
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with Mbh ^10^ Mq almost certainly reside in pseudobulges 
(e.g., Greene et al. 2008; Jiang et al. 2011). Between these 
two limiting cases (Mbh ~ 10^-10^ Mq), BH masses estimated 
from single-epoch spectroscopy presently cannot be known to 
better than a factor of ^ 2 without knowledge of the bulge type 
of the host galaxy. 
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TABLE 1 : Sample of Reverberation-mapped AGNs 


Name 

log ALa(5100A) 

FWHM(mean) 

log VP (mean) 

logMBH(RM) 

Bulge 


(erg s"^) 

(km s“^) 

(Tfo) 

(Me) 

Type 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

3C 120 

43.96+0.06 

1430+16 

0.291^3 

7.89l“Ti 

GB 

3C 390.3 

44.36+0.03 

13211+28 

9 40 + 0.02 

^•"*O_0.02 

9.24l“T 

GB 

Ark 120 

43.92+0.06 

6042+35 

1 KO+O.OS 
l.o^_q 03 

091 +0.14 
0-41_o.i6 

GB 


43.57+0.10 

6246+78 


s.iolg'.U 

GB 

Arp 151 

42.50+0.11 

3098+69 

0 . 181^6 

6.841“;}? 

GB 

Fairall 9 

43.92+0.05 

6000+66 

1 +0-03 

-L-01_o.o3 

0 cq+0.95 

0.00_o.i4 

GB 

Mrk 79 

43.57+0.07 

5056+85 

1 10-^0.04 
-L-J-O_0.04 

0 1 r + 0.16 

O-J-0_o.29 

GB 


43.67+0.07 

4760+31 

1 1 q+0.04 
-L-J-O_0.04 

8 . 021 “;}} 

GB 


43.60+0.07 

4766+71 

4.1 + _o.04 

7 oc:+0.21 

AoC)_o.i 2 

GB 

Mrk 202 

42.21+0.18 

1471+18 

U.OZ_Q JO 

5.981“;}} 

PB 

Mrk 509 

44.13+0.05 

3015+2 

i-03tg:gi 

0 1 q+0.10 

o.io_o.io 

GB 

Mrk 590 

43.07+0.11 

3729+426 

0 . 65 tVil 

7.651“;}} 

PB 


43.32+0.08 

2744+79 

0.5+Hi 

7.791“;}} 

GB 


43.59+0.06 

2500+43 

0.58^“:°t 

y yo+0.13 

1 • <d_o.i3 

GB 

Mrk 1310 

42.28+0.17 

2409+24 

U.lO_0 09 

6.501};}? 

GB 

NGC 3227 

42.36+0.03 

4445+134 

0.421^3 

7 40 + 0.27 
t •^'5-0.27 

PB 

NGC 3783 

42.55+0.18 

3770+68 

0.38l“:l“ 

7.151“;}} 

PB 

NGC 4253 

42.57+0.13 

1609+39 

-0.351H? 

5.981};}} 

PB 

NGC 4593 

42.79+0.18 

5143+16 

0 . 78 tVio 

6.831};}} 

PB 

NGC 4748 

42.55+0.13 

1947+66 

— 0 

'-'•-‘-^-0.08 

6.351};}} 

PB 

NGC 5548 

42.95+0.11 

12771+71 


8.091};}} 

GB 

NGC 6814 

42.08+0.29 

3323+7 

0 . 02 l“;« 

+ U/_o.ii 

PB 

NGC 7469 

43.36+0.10 

1722+30 

0 10+0.06 
'-».LO_o.o 6 

7 iq-I-0.13 
+ J-y_o 13 

PB 

PG 0003+199 

43.70+0.06 

1792+3 

n oc+O.OS 
U.OD_o.03 

( .+U_o.i2 

GB 


43.78+0.05 

1679+2 

n OO+0.03 

u.oo_o.o3 

7.161};}} 

GB 


43.68+0.06 

1273+3 

U.U41:_0 03 

7.461};}} 

GB 

PG 0026+129 

44.91+0.02 

2544+56 

1.301^2 

0 cc+O.lS 

0.00_o.i4 

CB 

PG 0052+251 

44.75+0.03 

5008+73 

1.801““2 

8.771};}} 

CB 

PG 0804+761 

44.85+0.02 

3053+38 

1 42+0.02 
i.^Z_ 0.02 

8.751};}} 

CB 

PG 0844+349 

44.18+0.07 

2694+58 

U.yO_0.04 

7.501};?} 

CB 

PG 0921+525 

43.62+0.04 

1543+5 

n 17 + 0.02 

j-1 _0.02 

* .UO_o.26 

CB 


43.69+0.04 

1658+3 

n 07+0.02 
u.z { _0.02 

i A'J_0.14 

CB 


43.47+0.05 

1600+39 

n 10+0.03 
'J-J^D_0.03 

' •-*-' 0 - 0.12 

CB 

PG 0953+414 

45.13+0.01 

3071+27 

1 KQ + O.Ol 

i. 0 o_ 0 . 0 l 

0 rq+0.12 
O.OO_0.i2 

CB 

PG 1211+143 

44.69+0.08 

2012+37 

0.98l“:°5 

8.02I};}} 

CB 

PG 1226+023 

45.90+0.02 

3509+36 


8.961};}} 

CB 

PG 1229+204 

43.64+0.06 

3828+54 

n 07+0.03 

A'.y < -0.03 

rj ry-i -j-O.ll 
* • ‘ -*-- 0.11 

PB 

PG 1307+085 

44.79+0.02 

5059+133 

-| OO+0.03 
i.OO_0.03 

8.891};}} 

CB 

PG 1351+695 

43.64+0.08 

5354+32 

1 07-I-0.04 
i _0.04 

7 Qc:+0.11 
+ oO_o.ii 

PB 

PG 1411+442 

44.50+0.02 

2801+43 

4.10_o 02 

0 90+0.16 
o.z,o_o.i3 

CB 

PG 1426+015 

44.57+0.02 

7113+160 

2 01+0 02 
^•'-'.L_ 0.02 

Q 07 + 0.12 
_o.i 2 

CB 

PG 1434+590 

43.73+0.05 

4711+49 

L.^u_0.03 

7 OQ+0.11 

I .oy_o.ii 

PB 


43.61+0.05 

5237+67 

1 qo+0.03 

L.^o_o.o3 

7.61I};}} 

PB 


43.61+0.05 

4767+72 

1 1 c:+0.03 
-L-J-'J_o.o3 

7.81I};}} 

PB 

PG 1613+658 

44.71+0.03 

9074+103 

^.^D _0 02 

0 c-i +0.21 

0 -Ol_o .21 

CB 

PG 1617+175 

44.33+0.02 

6641+190 

1 09+0.03 

L.o^_o.o3 

0 07 + 0.19 

o-o<- 0.12 

CB 

PG 1700+518 

45.53+0.03 

2252+85 

1 i:::9+0.04 

1 . 0 Z_o 04 

8.951};}} 

CB 

PG 2130+099 

44.14+0.03 

1781+5 

0.58l“;“2 

7.861};}} 

PB 

SBS 1116+583A 

42.07+0.28 

3668+186 

O.WtVit 

6.551};}} 

PB 


NOTE.— Col. (1): Name. Col. (2): Monochromatic luminosity at 5100 A for the AGN, corrected 
for host galaxy light and Galactic extinction. Gol. (3): FWHM of broad H/1, measured from the 
meau spectrum. Gol. (4): Virial product calculated usiug FWHM. Gol. (5): RM-based BH mass, 
calculated from the H/3 lag, line dispersion of broad H/3, and / = 6.3 for classical bulges and / = 3.5 
for pseudobulges. Col. ( 6 ): Bulge type: CB = classical bulge or elliptical; PB = pseudobulge. From 
compilation of Ho & Kim (2014), which also contains references to original data sources. 






